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We discuss a generalized Ising action containing nearest neighbour, next to nearest neighbour and plaquette 
terms that has been suggested as a potential string worldsheet discretization on cubic lattices by Savvidy and 
Wegner. This displays both first and second order transitions depending on the value of a "self-intersection" 
coupling as well as possessing a novel semi-global symmetry. 



1. INTRODUCTION 

Savvidy et al. jl| recently suggested a novel 
discrctizcd random surface theory, the so-called 
Gonihedric string, 

S= IJ2 \Xi~Xj\0M, (i) 

<ij> 

where the sum is over the edges of some trian- 
gulated surface, 0(otij) = \n — otij\^, C is some 
exponent, and ajj is the dihedral angle between 
neighbouring triangles with common link < ij >. 
This definition of the action was inspired by the 
geometrical notion the linear size of a surface, 
originally defined by Steiner. 

In equ.(jl]) it is the surface itself that is dis- 
cretized, rather than the space in which it is em- 
bedded. An alternative approach to discretiz- 
ing the linear size is to discretize the embed- 
ding space by restricting the allowed surfaces to 
the plaquettes of a (hyper)cubic lattice. This 
method was applied by Savvidy and Wegner |^- 
who rewrote the resulting theory as a gener- 
alized Ising model by using the geometrical spin 
cluster boundaries to define the surfaces. The 
energy of a surface on a cubic lattice is given ex- 
plicitly in the Savvidy- Wegner models by E = 
ri2 + Ann±, where n-i is the number of links where 
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two plaquettes meet at a right angle, 714 is the 
number of links where four plaquettes meet at 
right angles, and k is a free parameter which de- 
termines the relative weight of a self-intersection 
of the surface. In the limit k — > 00 the surfaces 
would be strongly self-avoiding, whereas the op- 
posite limit k — ► would be that of phantom sur- 
faces that could pass through themselves without 
hindrance. This should be contrasted with the 
standard 3D Ising model with nearest neighbour 
interactions where the surfaces are weighted by 
their areas. 

On a cubic lattice the generalized "Gonihedric" 
Ising Hamiltonian which reproduces the desired 
energy E = n-i + Ann^ contains nearest neigh- 
bour (< i,j >), next to nearest neighbour (<< 
i,j >>) and round a plaquette ([i,j,k,l]) terms 

<ij> 

+ H <J ^O a k(Tl- (2) 

[i,j,k,l] 

Such generalized Ising actions, or their equiva- 
lent surface formulations, have quite complicated 
phase structures for generic choices of the cou- 
plings |6[-^) . The particular ratio of couplings in 
equ.(|^), however, is non-generic and introduces a 
novel symmetry into the model - it is possible to 
flip any plane of spins at zero energy cost. 
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Although not a local gauge symmetry, the flip 
symmetry of the model is intermediate between 
this and a global symmetry. This symmetry poses 
something of a problem when carrying out simu- 
lations, as it means that a simple ferromagnetic 
order parameter 



(3) 



will be zero in general, because of the layered na- 
ture of the ground state. Even staggered magneti- 
zations would not do as the interlayer spacing can 
be arbitrary. It is possible, however, to force the 
model into the ferromagnetic ground state, which 
is equivalent to any of the layered ground states, 
with a suitable choice of boundary conditions on 
a finite lattice. Fixed boundary conditions would 
do the job by penalizing any flipped spin planes 
by a boundary term, but they carry a higher price 
in finite size effects than the customary periodic 
boundary conditions. A more elegant solution is 
to fix any two internal planes of spins in the lat- 
tice, whilst retaining the periodic boundary con- 
ditions. This has the desired effect of picking out 
the ferromagnetic ground state, whilst minimiz- 
ing any finite size effects. With fixed spin planes 
we can therefore still employ the simple order pa- 
rameter of equ.([|). 

2. ZERO-TEMPERATURE AND MEAN 
FIELD 

As the Gonihedric model is a special case of the 
general action considered in Q we can apply the 
methods used there for both the zero temperature 
phase diagram and mean field theory. For the 
zero temperature case this involves writing the 
full lattice Hamiltonian as a sum over individual 
cube Hamiltonians 



<h3> 
1 — K 



22 <y%OjO k ai 

i,j,k,l] 



(4) 



and observing that if the lattice can be tiled 
by a cube configuration minimizing the individ- 
ual h c then the ground state energy density is 
6{) = mini h c . 



This approach reveals that a layered ground 
state with parallel layers of flipped spins perpen- 
dicular to one of the lattice axes and arbitrary in- 
terlayer spacing is degenerate with the ferromag- 
netic ground state for all k. This might have been 
expected from the flip symmetry of the hamilto- 
nian itself. In addition, at n — an extra ground 
state corresponding to diagonal flipped planes ap- 
pears. 

In the mean field approximation the spins are in 
effect replaced by the average site magnetizations. 
The calculation of the mean field free energy is 
an elaboration of the zero temperature approach 
in which the energy is decomposed into a sum of 
individual cube terms. The next to nearest neigh- 
bour and plaquette interactions in the Gonihedric 
model give the total mean field free energy as the 
sum of elementary cube free energies <j){m c ), given 
by 



<t>{mc) 



K \ s 
9 2-^, 



<i,j>CC 



4 ^ 

<<i,j>>CC 

1 - K 



mirrijra^mi 
[i,j,k,l]cC 

Y^X^ 1 + m i)ln(l + m t ) 

iCC 

(1 - rrii)ln(l - rm)] 



(5) 



where mc is the set of the eight magnetizations 
of the elementary cube. This gives a set of eight 
mean-field equations 



d<j>(m c ) 
drrii 



= 



(i=i.. 



(6) 



one for each corner of the cube. Numerical iter- 
ation of these equations shows a single transition 
from a paramagnetic high temperature state to a 
layered, or the equivalent ferromagnetic, low tem- 
perature state. The [3 C determined in this fashion 
decreases quite sharply with k. 

3. MONTE-CARLO SIMULATIONS 

In order to see how the full theory tallies with 
the zero-temperature and mean-field results we 
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carried out Monte-carlo simulations |9j for various 
k values and lattices of size 10 3 , 12 3 , 15 3 , 18 3 , 20 3 
and 25 3 . Periodic boundary conditions were im- 
posed in the three directions and three internal 
perpendicular planes of spins fixed to be +1. A 
simple Metropolis update was used because of 
the difficulty in concocting a cluster algorithm for 
a Hamiltonian with such complicated interaction 
terms. The program was tested on the standard 
nearest neighbour Ising model and the some of the 
parameters used in the generalized Ising models 
of H to ensure it was working. We measured the 
usual thermodynamic quantities for the model: 
the energy E, specific heat C, (standard) mag- 
netization M, susceptibility x an d various cumu- 
lants. 

We can clearly see the second order nature of 
the peak in the specific heat in Fig.l 
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Figure 1. The specific Heat for k = 1 on lattices 
of various size. 



Although extracting the exponent a from the 
scaling of the peak is not particularly reliable 
because of the constant term that appears in 



both the direct and finite-size scaling forms an 
approach using the energy measurements gives 
(a — l)/v = —1.3(2). Binder's magnetization 
cumulant and the scaling of the susceptibility x 
which is shown in Fig. 2 
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Susceptibility 
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Figure 2. The susceptibility for n — 1 on lattices 
of various size. 



give an estimate of v = 1.2(1). One also finds 
= 1.79(4) from the FSS of x- All these expo- 
nents, rather remarkably given that the model is 
defined in three dimensions, are close to the On- 
sager values of the two-dimensional Ising model 
with nearest neighbour interactions. The ques- 
tion of whether this is simply a numerical coinci- 
dence or a true equivalence awaits an answer from 
further numerical work. The exponents and the 
transition temperature vary little, if at all, as k 
is increased. 

Simulations of the k = model llOll 



H = 



2 ^ 

[i,j,k,l] 



(7) 



4 



show that this is a special case, displaying a first 
order transition. The transition still appears first 
order at k — 0.1 but softens rapidly as k increases, 
so the crossover to the second order behaviour 
seen at n = 1 is quite sharp. 

4. CONCLUSIONS 

We have investigated a class of Ising-like mod- 
els suggested by Savvidy and Wegner as a lattice 
discretization of a particular random surface ac- 
tion that contains no bare surface tension term. 
The models have a novel flip symmetry that al- 
lows entire spin pnaes to be flipped at zero en- 
ergy cost. We have found that the phase diagram 
displays only a single transition to a low temper- 
ature layered state that, as a consequence of the 
flip symmetry, is equivalent to a ferromagnetic 
ground state. For non-zero k the transition is 
second order, whereas at, and probably close to, 
k = a first order transition is manifest. The 
presence of some degree of self- avoidance, in the 
form of a non-zero value for k, would thus appear 
to have an important influence on the universal- 
ity properties of the Savvidy- Wegner/gonihedric 
models. 

It is also possible to define open surface vari- 
ants of the models described here, which now in- 
corporate link spins as well as vertex spins, and 
to define models with similar properties in higher 
dimensions. All of these merit investigation both 
from the string theory and statistical mechanical 
point of view. 
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